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ABSTRACT 

A puzzhng feature of the C/ian(ira-detected quasar jets is that their X-ray 
emission decreases faster along the jet than their radio emission, resulting to 
an outward increasing radio to X-ray ratio. In some sources this behavior is so 
extreme that the radio emission peak is located clearly downstream of that of the 
X-rays. This is a rather unanticipated behavior given that the inverse Compton 
nature of the X-rays and the synchrotron radio emission are attributed to roughly 
the same electrons of the jet's non-thermal electron distribution. In this note we 
show that this morphological behavior can result from the gradual deceleration of 
a relativistic flow and that the offsets in peak emission at different wavelengths 
carry the imprint of this deceleration. This notion is consistent with another 
recent finding, namely that the jets feeding the terminal hot spots of powerful 
radio galaxies and quasars are still relativistic with Lorentz factors F ~ 2 — 3. 
The picture of the kinematics of powerful jets emerging from these considerations 
is that they remain relativistic as they gradually decelerate from Kpc scales to 
the hot spots, where, in a final collision with the intergalactic medium, they 
slow-down rapidly to the subrelativistic velocities of the hot spot advance speed. 

Subject headings: galaxies: active — quasars: general — radiation mechanisms: 
nonthermal — X-rays: galaxies 
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1. Introduction 

The superior angular resolution and sensitivity of Chandra has led to the discovery of 
X-ray emission from a number of quasar jets. Schwartz et al. (2000) were the first to note 
that the X-ray emission from the knots of the jet of the superluminal quasar PKS 0637-752 
(Chartas et al. 2001), at a projected distance ~ 100 Kpc from the quasar core, is part 
of a spectral component separate from its synchrotron radio-optical emission and it is too 
bright to be explained through Synchrotron - Self Compton (SSC) emission from electrons 
in energy equipartition with the jet magnetic field (note however that in the innermost knot 
of some sources - e.g. 3C 273 (Marshall et al. 2001), PKS 1136-165 (Sambruna et al. 2002) 
- a synchrotron X-ray contribution is possible). These properties are apparently common 
to other quasars jets, as indicated by the mounting observational evidence (Sambruna et al 
2001; Marshall et al. 2001; Jester et al. 2002; Sambruna et al. 2002; Siemiginowska et al. 
2002; Jorstad, Marscher, & McHardy 2003; Siemiginowska et al. 2003; Yuan et al. 2003; 
Cheung 2004, Sambruna et al. 2004). 

To account for the level of the observed X-ray emission, Tavecchio et al. (2000) and 
Celotti et al. (2001) proposed that this is due to External Compton (EC) scattering of the 
Cosmic Microwave Background (CMB) photons off relativistic electrons in the jet, provided 
that the jet flow is sufficiently relativistic (F ~ 10) to boost the CMB energy density in the 
flow frame (by F^) at the level needed to reproduce the observed X-ray flux. 

In all these sources the radio-to- X-ray ratio increases downstream along the jet, an 
unexpected behavior, given that the cooling length of the EC X-ray emitting electrons (7 ~ 
a few hundreds) is longer than that of the radio emitting ones (7 ~ a few thousands) and 
comparable or longer than the size of the jet, which would lead to a constant X-ray brightness 

as far out as the hot spots, contrary to observations. More surprisingly, in some jets (e.g. 
3C 273 in Sambruna et al. 2001 and Marshall et al. 2001; PKS 1136-135 and 1354+195 in 
Sambruna et al. 2002; PKS 1127-145 in Siemiginowska et al. 2002; 0827+243 in Jorstad, 
Marscher, & McHardy 2003) the X-ray and radio maps are anti-correlated, with the X-ray 
emission peaking closer to the core, gradually decreasing outward, and the radio emission 
increasing outward to peak practically at the very end of the radio jet. 

To explain the reduction of the X-ray fiux along the jet Tavecchio, Ghisellini & Celotti 
(2003) suggested that the X-ray emission originates in a collection of micro-knots undergoing 
adiabatic expansion sufficient to produce the desirable electron cooling. However, this would 
suppress also the radio emission, leading to practically indistinguishable radio and X-ray 
morphologies, contrary to their observed spatial anti-correlation. An elegant suggestion by 
Dermer & Atoyan (2002) that the X-rays are synchrotron emission from electrons cooling 
in the Klein-Nishina regime naturally produces shorter sizes in X-rays than in radio. How- 
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ever, it also produces larger optical than X-ray sizes, contrary to observations, and it seems 
therefore not to be applicable in this particular context. 

Our view is that, although X-ray producing electrons are present throughout the jet, 
the X-ray brightness decreases because the CMB photon energy density in the flow frame 
decreases along the jet as a result of a decelerating relativistic jet flow. The decrease in 
the flow Lorentz factor leads to a decrease in the comoving CMB photon energy density 
and hence to a decrease in the X-ray brightness along the jet. At the same time, the flow 
deceleration leads to a compression of the magnetic field; this results to an enhanced radio 
emission with distance, which gets thus displaced downstream of the EC X-ray emission. 

Based on the radio, optical, and X-ray jet maps, we argue in this note, that powerful 
extragalactic jets are relativistic and gradually decelerating. In §2 we formulate the syn- 
chrotron and EC emission process from a decelerating jet flow and present our results. In §3 
we discuss our findings and touch upon some open issues. 

2. The Model 

We parameterize the jet fiow assuming that its radius r scales as r{z) = rof{z), where 
f{z) is a function with f{zo) = 1 and tq is the jet radius at a fiducial point zq along the jet. 
Similarly, we assume that the fiow decelerates as T{z) = rQg{z), where g{z) is a function 
with g{zo) = 1 and Tq is the bulk Lorentz factor at zq. The evolution of the electron energy 
distribution (EED) along the jet is determined by the combination of adiabatic changes 
(losses for expansion and gains for compression), radiative losses, and particle acceleration. 
If the radiative losses of the electrons responsible for the X-ray and radio emission are 
neghgible (as may be the case for the (7 ~ a few hundreds) EC X-ray emitting electrons, 
or if we assume that some localized (e.g. shocks; Kirk ct al. 2000) or spatially distributed 
(e.g. turbulence; Manolakou, Anastasiadis, & Vlahos 1999) particle acceleration mechanism 
offsets the radiative losses, then the evolution of the EED can be treated as adiabatic. In this 
case, assuming that F ^ 1, the elementary volume dV scales as dV = dVof'^g, allowing for 
expansion or compression, while the adiabatic electron energy change rate with z is obtained 
from 7' = —7 (2/'// + g'/g) /3, where the prime denotes differentiation with respect to z. 
The solution of the above equation is 7(2;) = ^of~'^^^g~^^^- Assuming particle conservation, 
n(7, z)dVd'y — 71(70, Zo)dVod'yo, the comoving EED ^(7, z) can be written as 

n{j,z) = kr'r'^'^'^^'9-'^'^'^^', (1) 

where ^(7, zq) — k^~^ is the EED at zq. Following Georganopoulos, Kirk & Mastichiadis 
(2001), the EED in the local rest frame is ^(7, 2;)A^+*, where /^{z) is the familiar Doppler 
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factor A{z) — l/{r{z){l — (5{z) cos6)) and 9 is the observing angle. Using this, and taking 
into account the cosmological corrections, the EC flux per dz in the 5-function approximation 
(Coppi & Blandford 1990; an electron of Lorentz factor 7 upscatters seed photons of energy 
eo only to an energy 4eo7^/3) is 

^ karCU{l + Zy 2{s-l)/3 -is+2)/3A2+s f ^\ (2) 

dedz 2d^€o •' ^ \AeoJ ' ^ ^ 

where Z and d are the redshift and proper distance of the source, eo and e are respectively 
the energy of the CMB seed photons and the observed photons in units of mgC^ at Z = 0, 
and U is the CMB energy density a,t Z — 0. For 9 — 0, A{z) = 2T{z) and the z-dependence 
of the flux scales as oc /-2{5-i)/3^(4+2s)/3 

Assuming flux conservation, the magnetic fleld B{z) can be written as B{z) = Bq/ fg. 
The synchrotron flux per dz in the 5-function approximation (an electron of Lorentz factor 
7 produces synchrotron photons only at the critical synchrotron energy B^'^/Bcr) is then 



dF, 



1 R r(«+1)/2 

_ /gg^Cijer^o .-(7^-i)/6 -(5s+7)/6a(^+3)/2/ o \-{s-l)/2 /oN 

dedz" 127r(i2(l + Z)(^+i)/2^ ^ ^ ^'^^'■^ ' ^"^^ 

where Bcr — mlc^/eh — 4.414 10^^ G is the critical magnetic fleld. Note that while the 
synchrotron flux decreases with increasing redshift (the quantity d^(l + Z)^*"*"-*^)/^ increases), 
the EC flux (oc (1 + Z)'^/d'^) after an initial decrease up to Z Ri 1.9, remains practically 
constant for higher Z, as was flrst noted by Schwartz (2002). 

For ^ = the 2;-dependence of the flux scales as oc /~('''*~i)/6gr-(«-i)/3_ xhe condition 
that at ^ = the synchroton flux increases along the jet is that deceleration dominates over 
jet opening, leading to —g'/g < (7s — l)/'/2(s — 1)/. Taking the ratio of equation (3) to 
equation (2), the condition that Ls/Lec (or equivalently the radio to X-ray spectral index 
Urx) increases with z at ^ = is —g'/g < f j^j ^ and for a given jet opening profile / a 
milder deceleration g is needed to produce an increase in a^x than to separate the X-ray and 
radio peaks. This is in agreement with observations that show that the dechne of the radio 
to X-ray flux is a more widely observed trend that the X-ray - radio peak displacement. 

For the rest of this work we consider conical jets with / = (1 + az)/(l + a-Zo), where 
a = corresponds to a cylindrical flow, and deceleration proflles g = (z/zq)'", e > 0. The 
behavior of the EC and synchrotron flux proflles depends also on the observing angle. To 
demonstrate this, we plot in the upper (lower) panel of Figure 1 the EC (synchrotron) flux 
along a decelerating jet, for a range of angles. It can be seen that, while for small angles 
the EC flux peaks at the base of the jet, at larger angles the peak shifts downstream due to 
the fact that the emission from the faster base of the jet is beamed out of the line of sight 
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(for a given ^, A{z) peaks at r{z) ^ 1/sin^). The behavior of the synchrotron flux is also 
interesting. At angles 6^0 the peak of the synchrotron emission shifts downstream and for 
a given angle it is further downstream than the EC peak, in agreement with observations. 

It is interesting to see how much the above considerations change when radiative losses 
are important. The simplest case is that of a perfectly coUimated (i.e. a — 0) cyhndrical jet. 
This is also the case that requires the mildest deceleration for producing a given a^x increase 
or X-ray - radio separation. Assuming an EED slope s = 2, the EED remains a power law 
of the same slope following equation (1), with a maximum energy electron Lorentz factor 

e/3 

^ (z) - lmax{zo)Z^ , . 

"^"^"^^ ^ 1 + 7...(.o)[c^i(4^""^^/^ - 1) + cM'-''^' - 1)] ' 

where 7^0^(20) is the upper cutoff of the EED at Zq, = z/zq, Ci = 3zQABQ/87r{10e + 3), 
C2 = 4:ZqAU{1 + Z)'^Tl/{3 — 2e), and A = Aar/^nieC^TQ. The maximum EC energy eEc,max 
as a function of z is then eEc,max = ^o^"^ {z)j^^^{z) , while the maximum synchrotron energy 
is es,max — BozlA{z)j^ax{z) / ^cr- For a given observing energy and observing angle, one 
can find the radiating part of the jet and perform the optically thin radiative transfer by 
integrating expressions (2) and (3) along each line of sight for a two dimensional array of 
hues of sight that covers the source as projected on the observer's plane of the sky. 

The results of such calculations are shown in Figure 2, where it can be seen that the 
optical emission is confined at the base of the flow, due to the strong radiative losses of 
the high energy electrons, practically marking the site of strong particle acceleration. The 
morphology of the radio emission is angle dependent: at ^ = 6° it covers the entire extend 
of the jet, while at 6 = 12° the emission of its inner part is is dimmed because most of the 
radiation is beamed outside our line of sight. The X-ray emission at both angles peaks close 
to the base of the flow, due to the decrease of the comoving CMB photon energy density 
away from the flow base. Regarding the relative radio- X-ray morphology, while at = 12° 
the offset between the two images is very large, at = 6° the two images overlap, although 
the radio image clearly extends further downstream. Also, as the radio- X-ray spectral index 
plots show, in both cases arx increases downstream in agreement with observations. 

We compare now this picture with the one we would get from the same system if 
radiative losses were balanced by distributed reacceleration. (Figure 3). The X-ray images 
remain practically the same, because the radiative losses for the EC X-ray emitting electrons 
in this Z — 0.3, Fq = 6 jet are negligible. The radio image however changes significantly, 
particularly so for the ^ = 6° orientation. Now that the radio emitting electrons retain their 
energy the outer part of the jet is much brighter than the inner part, resulting to a clear 
separation between the radio and X-ray images. 



-6- 



3. Discussion and Conclusions 

We have proposed that the increase of the radio-to-X-ray flux ratio along the length 
of the jets of powerful quasars, as well as the occasional offset of the jet images in these 
wavelengths, are naturally accounted for in terms of relativistic flows that decelerate over 
the entire length of the jet. Despite this deceleration, the jets remain relativistic (F ~ a 
few) to their terminal hot spots (Georganopoulos & Kazanas 2003; hereafter GK03), within 
which they eventually attain sub-relativistic speeds. Our proposal provides, for the first time, 
a means for deducing the jet kinematics through simple models of their multiwavelength 
images; these can then be checked for consistency when coupled to the detailed hot spot 
emission, which as shown in earlier work (GK03), depends on the value of T at this location. 
We note here that a change of the Doppler factor can be also produced if the jet curves 
monotonically away from the line of sight without actually decelerating. However, in this 
case the CMB comoving photon energy density would remain constant along the jet, and 
the X-ray and radio emission would decrease along the jet in a similar manner, contrary to 
observations. 

The scenario we propose here finds further support from modeling jets with several 
sequential individual knots: Sambruna et al. (2001) noted the need for a gradual decrease 
of the Doppler factor and/or an increase of the magnetic field in order to reproduce the 
emission from the knots along the jet of 3C 273 with simple one zone models. These knot to 
knot variations can be naturally incorporated within the context of a decelerating coUimated 
fiow, as we propose. 

These same maps indicate also the need or not of distributed particle reacceleration along 
these jets, when the EC loss length scale of their radio emitting electrons (proportional to 
(1 + Z)^^r~^) is shorter than the observed length of the radio jet. This appears to be the 
case with sources like PKS 1127-145 aX Z = 1.187 (Siemiginowska et al. 2002) and possibly 
GB 1508+5714 at Z = 4.3 (Yuan et al. 2003; Siemiginowska et al. 2003; Cheung 2004), 
which show their peak radio emission displaced from that of the X-rays: in the absence 
of reacceleration, EC losses of the radio emitting electrons would produce radio jets shorter 
than the X-ray ones. Much about the reacceleration process at these jets can be inferred from 
their morphology at different frequencies. For example, the knotty optical jet morphology 
shows that reacceleration is not strong enough to offset the EC-dominated losses of the 
optically emitting synchrotron electrons (7 ~ 10^~^). 

The radiative efficiency of these jets is less than a few percent (e.g. Tavecchio et al. 
2000) and the energy lost in deceleration must be either used to heat up the matter in the jet, 
or must be transferred to material that is entrained by the jet. While the first option would 
result to an expansion of the jet, contrary to what is seen, entrainment from an external 
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medium would load the jet with baryonic mass while decelerating it. Entrainment would 
produce velocity gradients across the jet, and in this sense a faster spine and a slower sheath 
are to be expected. However, the observed X-ray and radio jet morphologies suggest that the 
dominant effect must be a deceleration along the jet. A fast spine that does not decelerate 
substantially and carries most of the jet power would produce a constant X-ray EC and 
radio synchrotron flux along the jet, in disagreement with observations. A consequence of 
entrainment would be that even if the jet did not start as a baryonic one, entrainment would 
gradually enrich it with baryons and eventually a fraction AF/F of its power, where AF is 
the decrease of the Lorentz factor, would be carried by the entrained baryonic matter. This 
in turn would increase the radio lobe equipartition energy content estimates derived under 
the assumption of a Icptonic composition. Our findings point to a picture where powerful 
relativistic jets decelerate, depositing most of their power in their surroundings in the form 
of kinetic energy, their observed radiation being only the tip of the iceberg. 

We want to thank the referee for his/her suggestions. G. M. thanks Rita Sambruna and 
Jessica Gambill for illuminating discussions on recent jet observations. 
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Fig. 1. — The EC (upper panel) and synchrotron (lower panel) flux normalized to its maxi- 
mum value for a range of observing angles along a jet that decelerates from Fq = 6 to F = 2 
with 6 = 1 from zq = 100 Kpc to Zmax = 300 kpc, while its radius doubles. The observing 
angles are ^ = 0° (short dash line), 9 = l/2Fo ^ 4.8° (long dash hue), ^ = I/Fq 9.6° (solid 
line), e = 2/Fo 19.5° (dot dash line). 
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Fig. 2. — The radio (top), optical (middle), X-ray (bottom) maps and the radio to X-ray 
spectral index a^x (below the X-ray map) for a decelerating flow observed under Q = '0° 
(left) and Q = 12° (right). The cylindrical flow, assumed to he at Z = 0.3, decelerates 
from To = 6 to r = 2 with e = 1 and Bq = 3. 10~^ as it propagates a distance of 200 
Kpc. The EED is a power-law of index s = 2 and 'jmax = 10^ at the base of the flow. 
This calculation includes radiative losses. Each map is normalized with red indicating the 
maximum luminosity per unit area. The projected jet lengths are ~ 21 Kpc for ^ = 6° and 
^ 42 Kpc for 9 = 12°. 
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Fig. 3. — The same as in figure (2) but with radiative losses balanced by distributed reac- 
celeration and without the optical map, which in this case would be identical to the radio 
map. 



